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ABSTRACT 


This study was an evaluation of a radiational scheme 
for inclusion in a numerical weather prediction model. 
This scheme uses empirical formulations for atmospheric 
absorptivities, scattering and reflectivity and cloud 
and earth surface reflectivity to compute the solar in- 
solation received and absorbed by the earth and key atmo- 
spheric layers. Longwave cooling by ae is also 
treated with empirically derived emissivities for water- 
vapor and carbon dioxide. This radiation model includes 
a two-layer cloud parameterization scheme and calculates 
the effect of an additional cloud layer on previous 
studies. 

In testing the validity of this model, use was made 
of Fleet Numerical Weather Central data in the form of 
vertical soundings on a spring day over the ocean. This 
data was used to calculate the radiation balance at the 
EHOpOpausce, tie caren s surtace and the key layers of 


the atmosphere. 





If. 


Pil. 


TABLE OF CONTENTS 


INTRODUCTION - =- - = ser r fr fr fF ke ke es ere 


DATA PREPARATION - - - = = = ef - = = feee- 


TEMPERATURE TERMS- - - - = -- 2+ eff ee 


MOISTURE TREATMENT IN THE SOUNDING - - - 


PRESSURE-SCALED ABSORBER MASSES- 


CLOUD PARAMETERIZATION = =~ = = = = += -- 


SOUR ORA DEAE LONG eee ee = = = ae ae ee ae ek 


A. 


B. 


PARTITION OF SOLAR INSOLATION= 


DISPOSITION OF F(S) INSOLATION 


i. 


2. 


3. 


DISPOSITION OF F(A) INSOLATION - 


1. 


2. 


Clear Sky Case - - - - - -=--e-e- 
Cloudy Sky Cases - - - - - ----- 


The Weighted F(S) Insolation - - - - 


The Clear Sky Case - - =- = = = = = = 
Overcast in both Cloud Layers- - - =- 


Disposition of F(A) Insolation with 
an Upper Overcast=- =~ = = \-=-feF-rte-- 


Disposition of F(A) Insolation with 
a Low OvercaSt - - = - - - - - - - 


Weighted F(A) Surface Insolation - - 
Atmospheric Layer Absorption of F(A) 
a. Overcast in both Cloud Layers- - 
De Upper Overcast Only- - - = - - = 
ce Lowere"G@vercast Only- =< —- - - = = 


d Composite Weighted Absorption 
by Layers- - - - ---+----- 


10 


Ls 


Be 


16 


iS 


21 


22 


22 


24 


Zo 


26 


27 


P5a | 


27 


30 


33 


33 


36 


37 


30) 


oi | 


eis) 


38 





IV. 


e. The Absorptivity 


(ABA) by Layers 


REFLECTION OF F(A) AND F(S) ENERGY - - - 


1. 


2. 


3. 


Albedo - - -\*-- 


Composite Absorptivity 


Check on Computations- 


STATISTICAL COMPUTATIONS 


1. 


2. 


Dee 


Clear Sky Cases- - 


Relationship between 
Cloudy and Clear Sky 


Statistical Relationship 


Albedo in the 
Cases - - -ef-o 


Ground Insolation in the 


Clear Sky Cases- =- 


Absorptivities in Air- 


General Conclusions-—- 


TERRESTRIAL RADIATLON= = = 


A. 


B. 


between 
Cloudy and 


THECRETICAL AND BEMPIRICAL BASIS= = = = <= 


x 
NED FLA te pWIth) CLOUDY SKIES = 


1. 


oa 


De 


NET FLUA F 


10 


* 
Net Flux Foo! with 


* 
Net Flux F - with 
in both Layers - - 


* e. ° 
Net Flux Foo with 


in the Upper Layer 


Net Flux F 
in the Low 


, wath 
Layer 


Oe * 
RO 


* 
a 
Composite Fio 


* 


6’ 


Clear Skies 


Overcast Clouds 


Overcast Clouds 


Only- 


Overcast Clouds 


Only- 


Calculations - - - - - 


With (eEOUDY SKIES. = = )=— = = 


* 
Wee or WITH CLOUDY SKIES = = = = = 


2 


BEE EEeCALTION, TO HEAT BAEANCE COMPUTATIONS 


Div lollLeaAL REsULTS— = 


te 


Downward Flux with Clear Skies - - - 


38 


39 


39 


42 


42 


42 


42 


44 


46 


47 


47 


48 


48 


51 


Syl 


[3 


23 


56 


10 


58 


63 


68 


68 


68 





Vi. 


2. Statistical Relationship between 
Sky and Clear Sky Cases- - - - - 


* * 
3. F. (CL) Relative to F . (O, O)- - 


MERIDIONAL CROSS-SECTIONAL DEPICTION OF 
RESULTS- - - = -=sefrerrrrer fr cfr ee 


CONCLUSIONS- - - = -efefrrrrrrtr ee 


INITIAL DISTRIBUTION LIST- 


FORM DD 1473 - - - ---+-+ef8e8 f+ 2 +282 20-2 


70 


71 


74 


82 


84 


87 


88 





1bdae 


LIST OF TABLES 


Ae Form of Original FNWC Grid Soundings- - - - 


b. Form of the Corresponding Radiative 
Soundings - - ---+-e+e8e 2222 22 -6e- 6 


Total Extraterrestrial Insolation (FADJ) at 
Each Grid Point of the Three Data Lines 
(ly minvl)- - - ---+e--+-+-ee-4-ee-e-e--8 - 


15 


40 





LIST OF FIGURES 
Five-Layer Radiative Sounding Model Used 
in this Study - - ----+-+-e--e----- - 


FNWC Polar Stereographic Grid and Meridians 
(1,2,3) Selected for Study- - - -----e-- 


Schematic Representation of F(A) Flux 
Disposition in the Clear Sky Case - - - - - = 


Schematic Representation of F(A) Flux 
Disposition in the Case of Two Overcast Layers 


Schematic Representation of F(A) Flux Dis- 
position with an Upper Overcast Layer Only- - 


Schematic Representation of F(A) Flux Dis- 
position with a Lower Overcast Layer Only - - 


* 
Terrestrial Net Flux F 


10 
a. Clear Sky Case- - - - - ---+---+-e-+ = 
Demereiwo OVerecastahayers (= S7= === i 


c. High Level Overcast - - - - = = ---e- 
d. Low Level Overcast- - - - - ----- = - 
* 

Terrestrial Net Flux Pe 

a. Clear Sky Case- - - ----+-+-+e---2e- 

Dee i vOmOVetedst lave row— 9 = =) = i= i= 

c. High Level Overcast - - - = = =e ff e- = 

d. Low Level Overcast- - - - ----+-fe-e-- 
* 

fhewres trial Net Flux F., 

a. Clear Sky Case- - - - - ----e-e--e- - 

b. Two Overcast Layers - ----ef#efrefefefe 2 

c. High Level Overcast - - - - - =~ ----- 


ar. WowWee OVe tea sta. = = = a eS eS Ss 


Jhal 


14 


28 


31 


34 


35 


D2 


54 


2 


a7 


59 


60 


61 


62 


64 


65 


66 


67 








Os 


itty 


Key to Mean Meridional Cross Sections- - - 


Mean Meridional Cross Sections 


ae Tropical Section - - - ------- = 


b. Higher Latitude Section- - - - - - = =- 





ACKNOWLEDGEMENTS 


The author wishes to express his appreciation to 
his wife for her support and for typing the rough draft 
of this thesis. 

Appreciation is also expressed to the author's 
thesis advisor, Professor F. L. Martin, for his sugges- 


tions, advice, guidance and support in this research. 





I. INTRODUCTION 


This study was an attempt to produce an empirical 
radiative parameterization for inclusion in the Fleet 
Numerical Weather Central (FNWC) primitiwe equation pre- 
aiection model. 

The radiative parameters described here include two 
cloud layers of fractional amounts CL(1) and CL(2) at 
each gridpoint instead of a single cloud layer as pre- 
viously used (Plante, 1973). The possibility of an addi- 
tional cloud layer tends to add substantial complexity to 
the radiation physics, but in general this paper repre- 
sents an extension of that described previously by Martin 
(1972) ¢ Some of the details of this radiation model are 
also based upon an unpublished manuscript (Martin, 1974), 
which in turn has some similarities to the multi-cloud 
models of Arakawa (1972) and Rodgers (1967). 

The vertical scale used in this study was based on 
the five sigma (0 = p/Tt) levels of the FNWC primitive 
equation model. The levels in this study are called k- 
eves (hoa 5 10) The k-levels are used for the 
purpose of data input aS shown in Figure l. Cloud cover 
is dealt with by means of a two-layer parameterization, 
with one high-level cloud layer from k = 4 to k = 6 and 


a low=Tevel cloud deck from k = 8 to k = 9, in amounts 


10 





Level- 


‘identifers O-level 
—————_ Onl 
Layer (0.2) i------- (AuqT,4e —------ 0.1 
| 2 o-0.2 
(2,4) 3------- Auqd de ------ 0.3 
| Ce aa o=04 
(4,6) 5-- --- - Sua, - —-—----—-- 0.5 
6—— 9206 
(6,8) a 4u,q.1,4e - -------07 | 
aa eae I ee 
(8 10) 9--------- 4u,q,T,4¢ - -----~-- 09 





o=10 


Figure 1. Five-layer radiative sounding used in this 
study. Levels are identified by their values on the k- 
scale, while layers are identified by their level boundary 
BiGieeou(eog., (a, LO) in parentheses,  Pressure-scaled 
water vapor and CO. mass increments AU and AC, respective- 
ly, as well as mixing ratio (gq) are introduced at the odd 
Mevels while temperature is introduced at all levels. 


Vea 





Gli) and €CL(2), respectively. In the computations for 
pressure-scaled water vapor and co. absorber masses a 
realistic method was used to extrapolate water vapor mix- 
ing ratio values above the level of k = 4, which is nor- 
mally the highest moisture level available from radiosonde 
data. The system for computation of the net flux of ter- 
restrial radiation makes use of new empirical absorptivi- 
ties of water vapor and co. (Sasamori, 1968), which are 
incorporated in the Martin (1974) Manuscript. 

To evaluate the radiative heat ee of this study, 
three meridians derived from FNWC analyses were used. 
These meridians are all from the Pacific Ocean area and 
most of the radiative soundings were in the Northern Hemi- 


sphere. 
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Il. DATA PREPARATION 


The data used in this study was in the form of three 
meridians of 25 vertical soundings over the Pacific Ocean. 
These meridians were located at 125°W, 170°W and 145°E as 
shown in Figure 2. 

The soundings on these data lines were taken from the 
Original FNWC analyses of temperature at the surface and 
at nine Since standard eentetaa Zo 100 millibars Lor 
OOOOGMT, 25 April 1973. The six moisture values were given 
by the vapor pressure at sea level, and by dew-point de- 
pression for levels from 925 mb to 400 mb. Table I(a) 
gives as an example the original sounding at gridpoint 
(1, 1) in the FNWC map. The data-processing methods 
applied to the temperature and mixing ratio data for the 


determination of the water vapor and CO, absorber masses, 


2 
and the cloud amounts within required k-level boundaries 
(Figure 1) are detailed in the following subsections. It 
should be noted that all soundings start at sea level, 
and as a reasonable approximation may be considered to 
have surface pressure 1 = 1000 mb. Hence the k-levels of 


Pguacaimenecame OOO, 900, S00, 700, 600, 500, 400, 300, 


200, and 100%mb, respectively. 


A. TEMPERATURE TERMS 
The temperatures were listed at each mandatory level 


Petween OC O07... L100 mb. The radiative Sounding 


Jue 





(63,1) 


((63-1), 1) 


(32,1) ine_< (32,32) 


» 


~ 
? 


jy h=10 





(1,1) 


Pugure 2.2), PNWC polar stereographic grid and meridians 
(1,2,3) selected for study. The longitudes A and hour 
angles of the sun relative to local noon are shown for 
each meridian. : 
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TABLE I(a). Form of original FNWC grid soundings. The 
humidity parameter between O95) M2., 400 -mb is dewpoint 
depression. 


Pressure (mb) T (2c) Humidity Parameters 
1000 28.5 27.7 mb 
925 20.0 a2 Se 
850 16.2 229°C 
700 7el S205 
500 -9.6 4.7°C 
400 -21.4 260°C 
300 -36.4 
250 -45.3 
200 - -55.5 
150 -66.9 
100 -80.3 


TABLE I(b). Form of the corresponding radiative sounding 
with temperature and mixing ratio listed at k-levels. 
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Pressure (mb) Tt 10K) Mixing Ratio 

1000 299.0 18.35 g(kg)> 
900 Zoe 16.50 
800 286.6 12.97 
700 280.3 7.84 
600 272.6 4.75 
500 263.6 2.97 
400 251.8 DIN 
300 236.8 0.73 
200 A gi 0.23 
100 192.9 0.04 
0 18229 0.00 





temperatures at certain k-levels were obtained by inter- 


polation. The temperature T was set equal to the listed 


10 


air temperature. At the top of the atmosphere, the tem- 
perature was assumed to be isothermal from 100 mb to p=0O, 
with the mixing ratio also approaching zero at p=0Q. 

The interpolation scheme used was the three-point 


Lagrangian formula of the form: 


Mo T) (CP, -P,) (PL mPL) IS ((PO-P,) (PamPo)) 
+ T, C(P -PQ) (PL -PS))/((P-P) (PP aPS)) 2a} 


+ 7, ((P,-P) (PL -PL))/((P,-P,) (P,-P,)) 


Here Po lies above P where k is the pressure level de- 


fen 


sired for T.~interpolation, and PirP. lie below level k. 


For example, with k = 9, Py and P. were taken as 925 mb 


oO! Ty and T. 
and P 


and 1000 mb, respectively. The parameters T 


are temperatures at the same levels as P Pp 


QO” i [esd yas 


respectively. 
B. MOISTURE TREATMENT IN THE SOUNDING 

The moisture parameters from the original sounding 
(Table I(a)) were converted into the mixing ratios at each 
of the original sounding levels. At k = 10 the surface 
vapor pressure, See! was used to calculate the mixing 


ratio by use of Equation (2-2): 


: | Dea 
nee 62197 e_, /1000 (2-2) 


where ec, Son mie 
air 


i 





— 
— 7 
7 - 
—_ 
= ae 
a 
=_ 
-_ = 





In calculating the mixing ratios at the remaining 
original data levels Equation (2-3) was used (after p. 63, 


Fleagle and Businger, 1963): 


(A- 2) ioctl a 
g(p) = 6 ey ee aie exp T Tp 1000 (2-3) 


Pp 


Here q(p) is the mixing ratio at level p. 

A = 21.656 

B = 5418.0°K 

T = temperature at level p (°K) 

T = dew point at level p (°K) 

T-T = dew point depression (°K) 

The q-values were then interpolated to the k-levels using 
an Penetton Similar to Equation (2-1). The results of 
this interpolation are shown in an example, Table I(b). 

Accurate humidity information is not generally avalil- 


able at or above 300 mb. To find the q-values at k = 3, 


2 and 1, an extrapolation formulation was used 


d 
q(p) ete (2-45 





This formula is similar to that developed by Smith (1966). 
However, instead of using the surface values of pressure 
and moisture in Equation(2-4) use was made of the pressure 
and moisture values at the 500 mb level. Equation (2-4) 
WacmeoivedmGer the parameter A using the “best-fit” between 


the variables y = log a and x = log ea . This was 


5 


7 





done by using the six values of q of each original sound- 
ings corresponding to p = O00) 925, 850, 700, 500, and 
400 of Table I(a). The "best-fit" values of A were then 
used in Equation (2-4) to extrapolate the q-values to 
p = 300, 200 and 100 mb. 

The usual formula for a best-fit parameter A relating 
Y; and x. samples is given in standard texts on statistics 


(e.g-, Crow et al, 1960) 





i=6 
yo, AER 
= 
A= = 
2 
x. 
ai 
i=l 


A further statistic computed in connection with each i\- 


profile was the correlation coefficient 





The correlation coefficient for each A-profile was general- 

ly in the range .95 to .99, which tended to enhance the 

credibility of the usefulness of the A-profile technique. 
The mean values, Ne together with the standard devia- 


tion, O(A), based upon the 25 soundings on each meridian 


are listed below: 


_ Merraian Ll Meridian 2 Meridian 3 
r 3-143 3-67 12 2.931 
G(X) 0.799 0,915 1.094 
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Cs PRESSURE=-SCALED ABSORBER MASSES 

With the mixing ratios computed at each odd k-level 
(Figure 1), it is possible to compute the pressure-scaled 
water vapor absorber masses in each layer. In Figure l, 
these water vapor layer-masses were denoted by AU, but 
here we develop notation which is more representative of 


the layers under consideration. For example, the pressure- 


scaled water vapor mass in layer (8, 10) is given by 


aie, se 18 
_ dg 9 0.72 7 
Oe) rm TORS GSU. ae) 


where the pressure ratio raised to the exponent 0.72 (after 
Moller and Raschke, 1964) has been employed. For the layer 
(6, 8), an analog to Equation (2-5) was employed with Io 
replaced by qd, and Po by Pos A similar procedure was fol- 
lowed in the topmost layers (0, 2),(2, 4) and (4, 6). 
Equation (2-6) affords an example for the topmost layer 
(O, 2). 
q, Ap p 
1 1 O22 

PANO aes) = g I[Toi3.25- abe 

In each of the five layers Ap is 200 mb. 


The algorithm for computing the integrated water vapor 


mass above the earth's surface is then as follows: 


ley = 020 

U (8) =U C10) +°U(S, 10) 

U(6) = U(8) + U(6, 8) 

U (4) = U(6) + U(4, 6) 2-7) 
URC 2) = U(4) Fr 2 4) 

U (1) = U(2) too WU (Oz. 2) } 

U (0) = U(2) + SUCCO 2) 


a, 





The principle for the computation of the carbon dioxide 
scaled mass is similar to that for water vapor. Here, how- 
ever, it is customary to state the pressure-scaled mass in 
terms of the normal temperature and pressure (N.T.P.) 
volume over a om? of the earth's surface. Thus 
4 Ap Po Oe 72 


g p_ (013.25) eer 





In (2-8), the fraction 3.14 x 1074 is the constant mixing 


Ap 





ratio of carbon dioxide by volume, while is the N.T.P. 
Ss 
The quantity Pp. is 


thickness of the layer AP = P 


One 6° 
the constant density in a standard homogeneous atmosphere 
for which se = 273.16°K everywhere and Dee 1013.25 7 mb. 


With minor modification, Equation (2-8) becomes 


Pp 





= -4 AP ) 0.72 
C(8, 10) = 3.14 x 10  (H) (yor37a5) “Fors. 25? (2) 
Bott 5 
H = 7 =— 7 99 5-10 cm. 


All other layer thicknesses of carbon dioxide are similar 


in form to Cie, 10), tor example: 


AP 7 O 


7 a4 S728 
Ae eae ranimin | ons D5. Ey) 


Pie winteqwaeuoan for sthe layer masses of C(&), C(6), C(4), 
Cil2Z) 7 2C (i) wana C(O) tollow exactly the procedure outlined 


EOnrowater «Vapor, 
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D. CLOUD PARAMETERIZATION 

In the calculations of the fractional cloud cover 
amounts in layers (4, 6) and (8, 9) the relative humidi- 
ties and thus the saturation vapor pressure at levels 
k = 5 and k = 9 ange be available. The saturation vapor 
Pressures at each level were calculated using Equations 
(2-ll(a)) and (2-11(b)), and the relative humidity using 


Equations (2-12(a)) and (2-12(b)): 


(A- =) 
_ .62197e 5 
q. (500) = aaa, (2-11 (a) ) 
(A- =~) 
_ .62197e 9 _ 
q_ (900) = 500 (2-11 (b) ) 
where A = 21.656 
B = 5418.0K 
_ q(500) - 
RH (5) “q_ (500) (2-12 (a) ) 
_ 9(900) 7 
RH (9) q_ (900) (2-12 (b) ) 


The relative humidity values were then used to calculate 
the fractional cloud amounts at each layer using the fol- 


lowing formulation, which is similar to Smagorinsky's (1960): 


CLT) ZOO R HAD) Oe 7 CZ 135(a.)) 


CL (2) $3559 (Kn (9) )— 92.0 (2-13 (b) ) 
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Tia. SOLAR RADIATION 


A. PARTITION OF SOLAR INSOLATION 

In this study, following Joseph (1966), a solar con- 
Stant of 2.00 ly aeua = was assumed. This is the value at 
the top of the atmosphere; however, it iS assumed subject 
to a four percent attenuation by oxygen and ozone above 
the tropopause. Thus the effective solar constant assumed 
at level k : 2 in this study is eZ i aie: 

The effective solar insolation at the tropopause was 


computed from 


F(2) = s[5-)~* cos Z@ (3-1) 
m 
where 
S = effective solar constant 
Cos Z = cosine of the zenith angle 
R/R_ = ratio of the actual earth-sun distance 


to the mean earth=-sun distance, 


The ratio R/R_, is listed as a function of the Julian date 
in the Smithsonian Meteorological Tables (List, 1958). 

In this study based upon data for 25 April at OOOOGMT a 
value of R/R_, = 1.00606 was used. The cosine of the 


zenith angle was calculated by 


Cos Z = Sind Sind + Cosd Cosd Cos h (3-2) 
where @ = Latitude 
6 = Solar declination for April 25th = 12°56! 
h = hour angle 


PA 





The value of the declination was taken from Table 169 
(List, 1958), and the value of Sind was calculated using 
two different formulas depending on the data-line used 
for the computations. These values of Sind along lines 
one, two, and three are calculated in terms of the FNWC 


map coordinates (I,J) as: 


973 2a 2(32-1)7 


Lines 1,3 Sind = 5 (3-3a) 
973.752 + 2(32-I1) 
. 973.152 = (32-1) 7 

Line 2 Sind = ——-—_ = is= 515) 


2 
973.752 + (32-1) 


Here I is the abscissa distance on the FNWC grid, Figure 2. 
The soundings were all taken at OOOOGMT; therefore, the 
value for the hour angle h, can easily be calculated. For 
example, in this experiment solar noon existed at the 180th 
meridian and the hour angles h = 55°, 10° and 35° occurred 
along lines one, two, and three, respectively. 

Following Joseph (1966) a computationally simple parti- 
tion of the solar insolation was used for this study. 
Joseph divided the insolation at level 2, F(2), into two 
parts, one of which was subject to water vapor absorption 


DUtEMOt t@nkRayler1gh scattering. This portion is given by 


F(A) = 0.349F (2); (3-4) 
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the remaining portion of the insolation, 
F(S) = 0.651F(2), (3-5) 


is subject only to Rayleigh scattering (but not to water 
vapor absorption) in a clear moist atmosphere. The 
presence of two cloud decks in this study introduces 
cloud-reflectivity into both the F(A) and F(S) solar 
energy portions. 

The F(S) energy is composed of all wavelengths d SS 
0.9 Um since here water-vapor absorption is negligible. 
The F(A) portion consists of those wavelengths X > 0.9 um 
where absorption by water-vapor and carbon dioxide bands 
in the near infra-red is the dominént attentuation process 


ine eclear aly « 


Bs DISPOSITION OF F(S) INSOLATION 
Joseph (1966) found that empirical values of Rayleigh 
clear-sky albedo, (after Coulson, 1959), were best-fit 


uSing the least squares technique by the following formula 


a(R) = .085 + .25074 [log,, — Sec Z)] (eG) 


0 


Il 


Pp LGi3.25 


0 


Here T/Po = 1 in view of the fact that mean sea level 
occurs near 1000 mb. The value Sec Z is equal to (Cos ayo 
where Cos Z is calculated from (3-2). 

Another reflective-type parameter which must be con- 


sidered in this study is the surface albedo a(G). Over 
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the land areas of the Northern Hemisphere a constant 
average surface albedo of 0.14 is utilized. However, the 
value over the oceans depends on cos Z. Since in this 
study all of the soundings were over the ocean, the 


equation 
a(G) = max {.06, .06 + .54 (.7 - cos Z)} (3-7) 


Was used to calculate 0(G), (after Gates et al, (1971)). 

Since a two-layer cloud model was used, four separate 
cloud cases a discussed below. 

OF Clear Sky Case 

In this case the F(S) insolation is subject only 

to Rayleigh sky reflectivity, a(R), and the earth-surface 
reflectivity, @(G), of Equations (3-6) and (3-7). Con- 
Sidering the possibility of multiple reflections between 
the earth and atmosphere, each of which diverts downward 
the fraction Q@(R) of the earth's surface reflectance Q(G), 
the clear sky insolation at the earth's surface after 


sar aeereee is given by 
TS10(0,0) = F(S) [l-a(R))] [1-a(G)]/[1l-a(R)a(e)] (3-8) 


In the next three cases, where clouds were present, 
the cloud reflectivity values suggested by C. D. Rodgers 
(1967) were used. Thus for the higher level clouds (between 
k = 4 and 6), R(1l)= 0.54 and for the lower level clouds (be- 


tween k = 8 and 9), R(2) is 0.66. 


PR, 





Cee lloudy oy Cases 
We consider first the fraction of the grid area 
having clouds at both levels. The F(S) insolation ab- 


sorbed by the ground at such a grid point is given by 


Pome = oP (SILC laR(1)) (1l=R(2)) (1l-a(G))] (3~9) 
Pea CCR (Ri ee R2) (CG) +R i) a (G) ) 


2 R (Re een — 


(after Arakawa,1972). The special symbols (1,1) denote 
Clie ="Clhi2) = 1.0, which will temporarily be assumed, 
The other cloud=-cover variations to be considered 
are those of CL(1) = 1.0, CL(2) = 0.0 and CL(1) = O.O, 
CL(2) = 1.0, respectively. In the former case, R(2) = 


0.0 in (3-9) so that 
TS$10(1,90) = F(S) (1l=-R(1)) (1-a(G))/(1-R(1)a(G)) (3-10) 


In the latter case with clouds at the lower layer only it 


follows that 
oO (Or ets) ClaR(2)) (1-0 (6G))7 (1-Rt2)0(G)) =) 


Equations (3-8), (3-9), (3-10), and (3-11) were used to 
Calculate, the absorbed F(S) insolation at the earth's sur- 
face after only reflective-type processes are considered 
to take place. However, in each case the sky was con- 
Sidered to be either overcast or clear and an appropriate 
grid-area weighting scheme was necesSary to correct to the 


acrid cloud Gonditions CL(1), CL(2). 
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Se The Weighted F(S) Insolation 


The cloud combinations denoted by (0,0), (1,1), 
(1,0), and (0,1), respectively, actually occur with the 
areal coverage (or weights) at a grid point given by the 


following functions: 


W(O,0) = (1-CL(1)) (1-CL(2)) 

W(1,1) = (€L(1)) (cL(2)) | (3-12) 
W(1,0) = (CL(1)) (1-CL(2)) 

W(0O,1) = (1-CL(1)) (CL(2)) 


ac. e@ach Grid point. 

When these weightS were applied multiplicatively 
in sequence with the right sides of Equations (3-8), (3-9), 
(3-10), and (3-11), then the composite surface-absorbed 


insolation, IS10 resulted. 


£51 0= 1S70(0,0)W(0,0)+ I1810(1,1])W(1,1) 


ee od Ode, Oo )evrcls Ot FSlOo( Ot) WO, -1) (3-13) 


Note that by subtraction of IS10 from F(S) we also 
arrive atuthateportton of the F(S) iansolation reflected to 


Space. 


Cc. DISPOSITION OF F(A) INSOLATION 
In the following sections, the fraction of the insola- 
tion which is absorbed by atmospheric water-vapor and 
earbon Gioxzde is discussed. 
ae ine clear sky Case 
In this connection the Manabe-Moller absorptivity 


Pune t 10 nm 


ate pek) = 2271 (U(2,2k) Sec pares (3-14) 
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Figure 3. Schematic representation of F(A) insolation dis- 
POSIti1On 2n thercilear sky case. 
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was used. The form of Equation (3-14) indicates that the 
absorptivity a is applied over the pressure-scaled water- 
vapor mass between levels 2 and 2k (see Figure 1) along 
the zenith slant-path Ange Ze. Therefore, assuming negli- 
gible water-vapor mass above level 2, the absorbed solar 


insolation in the layer (2,6) is 


Al(2,;6) -)O0.2.271F (A) ([U(2,6) Sec as acacia (3=15)) 


Similarly, the absorbed insolation in the layer (2,10) is 


AZ, 10))="0-.271F (A) [(U(2,10)Sec Zee (3-16) 


and the absorbed insolation in the layer (6,10) is given 


(see Figure 3) by 
A(6,10) = A(2,10) - A(2,6) =) 


The direct transmission of F(A) insolation to the 
earth's surface was found by subtracting A(2,10) from F(A). 
This direct transmission is partially reflected by the 
earth's surface (the value computed for a(G) in Equation 
(3-7) is applicable here) and the remainder, absorbed at 


the earth's surface, 1S given by 


IA10(0,0)=F(A){1-.271[U(2,10)Sec 77 ial 


(1-a(G)) 
(3-13) 
In the following subsections dealing with partially 


cloudy layers, it is necessary to adopt representative 


cloud reflectivities. For this purpose, values suggested 


Ze 





Dy C. D. 


Rodgers 


(1967) were used. 


These reflectivity 


values for the F(A) wavelengths differ from those listed 


by Rodgers for the F(S) 


Proceqaure 


wavelengths. 


We followed the 


of considering overcast cloud conditions when- 


ever clouds occur and subsequently applying the weight- 


ing factors as was done in Section III(B). 


2. Overcast in Both Cloud Layers 
Here, with clouds at both levels the cloud reflec- 
In 


ivity values RA(1) = .46 and RA(2) = .50 were used. 


addition Rodgers suggests cloud absorptivities of A(l) = 


~20 and A(2) 


set of formulas, 


CEC. ,; 


-30 (see Figure 4). 


With the following 
the downward directed insolations FA4}¥, 


were calculated at each required level in order to 


arrive at a value for the insolation abosrbed by the 


atmospheric layers, 


A(2, 
FA4Y 
a2, 


No 


4)' 


Ne 
=) 


6) 


8) 


9) 


F(A) .271(U(2,4)See 2]°°?? 


F(A) (1-.271(U(2,4)Sec 2)°°797) 


RUIRAC Tyee pune Avisee Zz) °° 


Az 4) eto) 


FA4+4A (1) 


FA4Y (1-RA(1)-A(1)) 


FA6+(TD(6, 8))[1-RA(1)RA(2) (TD(6, 


FA6+ - FAS8Y 
FASYA (2) 


FA8+(1-RA(2)-A(2)) 


Baer D (oe lO) (| l=o(G)RA(2) (TD (9, 10))-) 


FAQ9+ - FA1OY 
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and finally by the earth's surface. 


(3-19a) 
(3-19b) 
(3-19c) 
(3-19d) 
(3-19e) 
(3-19£) 
yon 
(3-19q)) 
(3-19h) 
(3-191) 


(Sa J) 
Sl. 


(319g) 
(3-192) 
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Pgivec 4 EmGechematic representation of F(A) insolation dis- 
position in the case of two overcast layers. 


eye 





where 


TD(6, 8) =1- SOUL (Gsen Ble (3=20) 


Pee TTUuoeTONS 7a) 


TD(9, 10) 


In the preceding equations A(2, 4), A(4, 6), 
A(6, 8), A(8,9), and A(9,10) represent the insolation ab- 
sorbed in the indicated layers. A superscript prime fol- 
lowing A(2,4) indicates absorption (by use of Equation 3-14) 
of the downward directed beam, while the double prime indi- 
cates absorption of the beam specularly reflected upward 
EZom the cloud top. The functions of Equation (3-20) repre- 
sent transmissivity after diffuse passage through a cloud 
or sub-cloud layer, which causes the effective replacement 
St sec 2 in BQWatlon- (3-14) by Sec Z = 5/3 (Katayama, 1966). 
The symbols FA6Y, FA8Y, FA9¥, and FALOY indicate the down- 
ward streams of insolation passing through the indicated 
levels. Note that the effect of multiple reflections be- 
tween clouds is incorporated in the denominator of FA8\V, 
while that between the lower cloud and the ground has a 
corresponding term in the denominator of FA1O+, involving 
the reflectivities a(G) and RA(2). The effect of the 
multiple reflectances as modeled here is to enhance the 
insolation at the lowermost reflecting boundary of the two 
under consideration. 

BFinawty, note that the flux absorbed by the earth's 
surface under the conditions of an overcast layer at both 


levels is given by 


PAO Cl ) p= sPALOC 1 =2(G) ) (3-21) 


a2 


ee 





3. Disposition of F(A) Insolation with an Upper 


Overcast 

In this case with only an upper overcast layer, 
the values of RA({1) and A{l) used are the same as those 
of the preceding section. However, a somewhat simpler 


set of disposition equations result: 


Ae = E272 1U(2, 4) sec Zi 
eu PMSeez +s > 77 RA(1) } (3-22a) 
Cea CAA eneacu ee 4)isec ea): >) te. (3-295) 
BAG = mtAy ie eT (uo ea iseens) °° TLISRA(1) =A (1)1 
(S=22¢)} 
Olea 7 (oto onsy3) 2° 4 
*{l-a0(G)RA(1)TD(6, 10)7}7> (3-224) 
A(6, 10)= FA6+ ~ FA10+ (3-22e) 


Here the variables are defined similarly to those of the 
preceding section. (Note that A(2, 4)' and A (2, 4)" have 
been combined in Equation (3-22) to give A(2, 4) and the 
Sermnrtron, of TD(6, 10) is similar in principle to that 

Oty TD (6,8) );- The absorbed insolation at the ground is 


given by 


TA10(1, 0) = FAI1OY (1 - a(G)) (3-23) 


(see Figure 5). 
4. Disposition of F(A) Insolation with a Low Overcast 
The set of equations which describe the behavior 
of the insolation with a low overcast deck are given 


below (see Figure 6). 


a3 
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Figure 5. Schematic representation of F(A) insolation dis- 
POSsttiOn Withean Wpper overcast layer only. 
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FLgure 6, Schematic representation of F(A) insolation dis- 
position with a lower overcast layer only. 
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ee 





0.303 


PuGe hoje =. Z271F (A) [U(2, 6$)Sec Z] (3-24a) 
A(6, 8)' = .271F(A)[U(2, 8)Sec eee 
0.303 
-(U(2, 6)Sec Z) ] (3-24b) 
FA8Y = F(A) [1-.271(U(2, 8)Sec ay Comer (3-24c) 
FA9Y = FA8Y[1-RA(2)-A(2) ] (3-24d) 
FASt = FA8YRA(2) (3-24e) 
A(6, 8)" = FA8Y[.271(U(6, 8)Sec ae (3-24f) 
A(6, 8) = A(6, 8)' + A(6, 8)" (3-24g) 
A(2, 6)" = .~271FA8t[(U(2, 8)Sec py gee 
0.303 
-(U(6, 8)Sec 2Z) “| (3-24h) 
Pace op = eae, 6)": + AC2, 46)" (3-241) 
A(8, 9) = FA8Y A(2) (3-243) 
FA1LOY = FAQOYTD(9, 10) [1-a(G)RA(2) (TD(9, GY oe Pals 
(3-24k) 
A(9, 10)= FA9Y¥ - FA1OY (3-242) 


These variables also have Similar definitions to those in 
Section III-C-2. As before the absorption by the earth's 


surface 1s given by 


IA10(0, 1) = FA1OY (1 - a(G)) (3-25) 


5. Weighted F(A) Surface Insolation 
As previously described in the section on the 
F(S) insolation, the values of IA10 calculated for the 
Gases (05,70), 0(1, 1), (1, 9) and (0, 1) must be appro- 
priately weighted by the values described in Equation (3-12) 
in order to arrive at the composite F(A) surface-absorbed 


insolation. 
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TALO = ITAL0O(0, O)W(0, O)+ IALO(L, L)WCS, 1) 


+ IA1O(1, O)W(1, O)+ IAI1O(O, 1)W(O, 1) (3=26) 


GC Atmospheric Layer Absorption of F(A) 


In each of the four cases discussed, the F(A) in- 
solation was absorbed as it passed through the key atmos- 
.pheric layers (2, 6) and (6, 10). In the clear-sky case 
these absorptions are given by Equations (3-15) and 
(3-17). However, in the calculations of the cloudy-sky 
cases the various contributions have merely been added 
together to give the total insolation absorbed in the 
layer (2, 6) and then in the layer (6, 10), even though 
cloudiness affects the radiation physics differently within 
the cloud and subcloud layers as Figures 3, 4, 5, and 6 
indicate. 

are Overcast in Both Cloud Layers 

In this case the insolation absorbed in the 
two layers was calculated with values from the appropriate 


equations from (3-19a, ..., 3-192) 


A(2, 6), “ac A(2, 4)' + A(2, 4)" + A(4, 6) 


A(6, ay: 


I 
Me 
en 
0 


1 + A(8, 9) + C9 3 10} (3-27) 


Where the Sulposeripe (1, 1) on A(6, 10) indicates the cloud 
layer combinations which will subsequently be weighted. 
Ds Upper Overcast Only 
For the high level overcast only the insola- 


tion absorbed was calculated using the various equations 


aa] 





iS =2 2a oe ge See Ze). 


» 
oe, 
NO 
~ 
.e)) 
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A(2, 4) + A(4, 6) 


FAGY - FA1OY (3-28) 


hed 

oe, 

ey) 

= 

j_ 
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Cs Lower Overcast Only 
In the low level overcast case the value of 
the insolation absorbed in the layers (2, 6) and (6, 10) 


was calculated using equations of (3-24a,..., 3-242) 


A(2, 6) = & (2,26) > + AZ; 6)” (3-29) 


od 
om, 
ey) 
~ 
© 
© 
~ 
It 
> 
ey) 


8)" + A(6, 8)" + A(8, 9) + A(I, 10) 


d. Composite Weighted Absorptiom by Layers 
In this case, as in the previous subsection 
on weighted F(A) surface insolation, the values computed 
in each of the four cases must be appropriately weighted 
by the functions described in Equation (3-12) in order to 
arrive at a composite value of the F(A) insolation absorbed 


in the atmospheric layers (2, 6) and (6, 10). 


A(2, 6)=A(2, 6)9 9 W(O, O)+A(2, 6), , W(H, 1) 

+A(2, aL ae W(1l, O)+A(2, ood W(O, 1) (3-30) 
A(6, 10)=A(6, 10), ,W(O, O)+A(6, 10), ,W(1, 1) 

#A(6, 10), QW(1, 0)+A(6, 10) ,W(O, 1) (3-31) 


eo The Absorptivity (ABA) by Layers 
In this study it was convenient to deal with 


GiomiractrOnateabSorptivity, Or Simply a@®Ssorptivity, as 
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well as the actual insolation energies absorbed in each 
layer, given by Equations (3-32 and 3-33) and by Equation 
(3-26) at the ground. In the calculation of the absorp- 
tivity, the total incoming insolation at the top of the 
atmosphere must be known. This value was calculated us- 


ing Equation (3-32) 
FADJ = 2.00(R/R_)* Cos Z (3=32) 


Table II lists the instantaneous values of FADJ at each 
gridpoint along the three data lines. The calculation of 
the full-tropospheric absorptivity then follows from 

A(2, 6) + A(6,- 10) 


ABA = ADI (3-33) 


which is the ratio of the insolation absorbed in the tropo- 
Sphere to the extraterrestrial insolation at the same time 


for each gridpoint. 


D. REFLECTION OF F(A) AND F(S) ENERGY 
In the calculations of the albedo and the fractional 
absorptivity of the ground both the F(A) and F(S) frac- 
tions of the solar insolation had to be considered. 
1. Albedo 
That portion of the F(A) insolation which was re- 


fiilected back to space 1s given by 


REFA = [F(A) - A(2, 6) - A(6, 10) - IA10] (3-34) 
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TABLE II. Total extraterrestrial insolation (FADJ) at 
each grid-point of the three data lines (ly min7l). 


Line 1 Line 2 Line 3 
~899 1.899 1.346 
923 aS a ees 1.374 
947 1.923 1.403 
2970 1.933 1.430 
~995 1.939 1.458 

1.018 1.946 1.485 

e042 1.947 ia § 0) 

1.064 1.947 1.534 

1.085 1.942 1,553 

1.106 1.934 oy Ae, 

Le Z5 Nees ees 1.597 

US ars 1 eNO) F; 1.613 

1.158 1.888 1.626 

is yaa 1.862 Ls634 

eS 0 1.834 1.638 

is bS7 1.798 1,637 

1.190 e757 1.629 

ke} 710 1,616 

1.181 1.658 eos 

1.168 1.598 1.565 

1.149 1.534 1.526 

lle 2 1.462 1.478 

1.088 1.384 1.419 

1.048 es Ol 1.350 
~998 1.210 i270 
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where the calculation equations of A(2, @), A(6, 10), Beye 
IA10O have been described in Equations (3-30), (3-31), and 
(3-26), respectively. 

The F(S) ie olacion reflected was simpler to calcu- 
late because the F(S) wavelengths do not undergo atmospheric 
absorption in the passage through the layers. The value 


of the reflected F(S) insolation is givem by 
REFS = F(S) =-/1S1L0 (3-35) 


where IS10 was calculated by Equation (3-13). 
These values of REFA and REFS were then combined 
to give a value for REF, which is the total reflected energy 


passing outward through the tropopause (k=2) 
REF = REFS + REFA (3-36) 


The albedo is defined as the ratzo of the total re- 
flected energy to the total incoming energy. This. ratio 


was calculated by 


REF 
= 32357 
ALB J ( ) 


where FADJ is defined in Equation (3-32). The values of 
REFA and REFS were computed at each grid-point. Taeene— 
sulting REF values are shown in cross-se:ction form in 
Figures ll(a) and l11(b). Also, values o€ the albedo were 
computed. The global mean albedo over the 75 soundings 
turned out to be .412, based upon the cl@®ud-cover para- 


Mererizatcton for ©Cl()), CL(2). 
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2. Composite Absorptivity (ABG) by Exarth 


The total wnsolation absorbed ate the earth was 
the sum of the weighted values of the F(S) and the F(A) 


portions of the incoming solar insolation. 
QO(E) = IA10 + 1S10 (3-38) 


Here Q(E) is the insolation aincorbed by tthe earth's sur- 
face at each grid-point sounding. 

Using Q(E), a fractional absorptivity for the 
earth's surface was calculated. The equation used for 


the calculation was 


OCB) z: 
ABG = —~TaDg (3-39) 


where the calculations of Q(E) and FADJ have been previously 
discussed. 

3. Check on Computations 

The actual computational process used in this model 

was checked against the condition that the summation of the 
fractional values ALB, ABA, and ABG at each grid-point 
equals .96. This must be expected since we allowed for 
attenuation of the solar insolation in the amount of four 
DeTeent as, it passed through the stratosphere. This check 


computation was satisfied at each grid-point. 


E. SEAT ror realk, COMPUTATIONS 


ie Clear Sky Cases 


In Order to werify the preceding calculations, 


statistical computations of several key values were obtained 
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uSing linear regression techniques from the BMD set of 
statistical programs (Dixon, 1973). The predictand of 
the first computation was ALB (0,0,Z), which is the 
albedo in the clear sky case. The predictor was the se- 
cant of the zenith angle (Sec Z). The best-fit equation 


which resulted was 
AGB CO, Opa is =9(—0.02726)7 0.71387 Sects) (3-40) 


with the regression coefficient R= .979. Equation (3-40) 
summarizes how the radiative model linearizes the re of 
Sec Z in the clear sky case. The high R-value confirms 

the known result that Sec Z is a good predictor of the 
albedo (with clear skies), which is to be expected over 

an oceanic area, The mean albedo in the clear sky caSe, 
which corresponds to the average secant of the zenith 


angle, is as follows 
AGE® (0.0) = 5.235 Sec Z = 1.428 


Two other regressions that were performed using 
the linear regression technique were ABA (0,0,M) and 
ABG (0,0,M) against the water-vapor mass path length M, 


where M is’ given by 


{2 2 


M= (U Sec rie (logy, (US Sec z)7 ] 


This parameterization of the water vapor mass is essential- 


ly that of Hanson (1971), who used a similar absorber-mass 
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predictor in his empirical specification of ABA for both 
clear and partly cloudy skies. The resuklting equations 
obtained with clear skies in this study for ABA (0,0,M) 


and ABG (0,0,M) were 


/ [2 


-096 + .031(U Sec z)t J 


ABA (0,0,M) * (log, 9 (U Sec zy 


R= .984 (3-41) 


{2 


[2 1 
[log, ,(U Sec. 4) ] 


-800 - .035(U Sec zr 


ABG (0,0,M) 


R= .956 (3-42) 


The means and standard deviations of ABA(0,0,M), ABG(0,0,M), 
and M, based on all of the original 75 smundings, are 


given below 


ABA(0,0,M) = .137 G = 92029 
ABG(0,0,M) = .687 Oo = .034 
M = 1.342(gm ene ee 6 =m 054a 


2. Relationship between Albedo in the Cloudy and 
Clear Sky Cases . 


As discussed previously, it was possible to compute 
the value of ALB in the clear sky case and also in the com- 
posite cloudy sky case. In order to study the relationship 
between these two cases a ratio ALB[CL(1LJ,CL(2))]/ALB(0,0) 
for the same zenith angle at grid point (01,3) was formed. 
ipissvalwe was used as the predictand ane the total opaque 
cloud cover (CL), specified by the model was used as the 
predictor. Here CL is calculated using!@@uation (3-43), 


which approximates the effective cloud-cover by the two 
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layers of clouds: 
CL = CL(1) + CL(2) - CL(L) (CL(2) (3-43) 
The regression equation which resulted was 


ALB([CL(1),CL(2)] = ALB(0O,0) (1 + 3.931CL —- .8782CL“) (3-44) 


R = 0.956 


The means and standard deviations of the values in Equa-~ 


tion (3-44) are 


ALB[CL(1),CL(2)] = .412 o = .43]1 
ALB(0,0) = .135 o = .141 
cL = .680 GO = .742 


Note that the mean cloudy sky albedo generated by 
the model is 3.211 times that of the meam clear sky albedo. 
The high value of the correlation coefficient (R = .956) 
indicates reliability of the indicated dependence upon CL. 
This global mean albedo of .412 is somewhat Cee than 
that recently reported [ALB = .300] by Raschke et al (1973). 
However, the mean value reported here does not have an 
input from the polar latitudes. Raschke, et al suggest 
that perhaps the reduced cloud amounts (CL) recently re- 
ported in low latitudes would account for the reduction 
compared to that computed here by the parameterization 


Equation (2-13), after Smagorinsky. 
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3 Statistical Relationship between Ground Insolation 
im ene Cloudy and Clear Sky Cases 


In this section as in the preceding one a ratio 
Y = ABG(CL)/ABG(0, 0) was formed for each sounding case 
to determine the best-fit regression formula using CL as 


the independent variable. The formula which resulted was 
Y = ABG(CL)/ABG(0,0) = (1 — .66CL), R= .999 (3-45) 


The ratio of ABG(CL)/ABG(0, 0) iS essentially the 
corresponding ratio of the transmissions received by a 
pyrheliometer. Quinn (1970) has summarized the empirical 
solar transmissivity formulas by many investigators and 
shows that some of the most reliable estimates are very 
nearly as formulated by Equation (3-45). For example, 


Budyko (1956) for the tropics lists 
ABG(CL)/ABG(0, 0) = 1 — .655CL 

and Kimball (1928) lists the corresponding ratio as 
ABG(CL)/ABG(0, 0) = 1 - .71C 


The mean statistics over the sample soundings gave 


ABGUCE) = 2379 
ABG(0, 0) = .688 
CL = .680 
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an RDSOrPelvcles Tee glamey Gr « 


Here the computed values of the ratio ABA(CL)/ 
ABA(0O, 0) were investigated by least squares with the 


result 
ABA (CL) /ABA(0O, 0) = 1 + .37863CL, R= .781 (3-46) 


with mean statistics for the 75-case population 


ABA (CL) = .169 
ABA(O, 0) = .137 
CL = .680 


This indicates that the atmospheric absorptivity is some- 
what enhanced with the cloud-model solar dispositions in- 
troduced. An increaSing solar absorptivity with increasing 
cloud cover is in general agreement with Plante's results 
(1972) in a one-layer cloud model. 
5. General Conclusions 

The statistical relationships developed from the 
radiation model show consistency with results reported by 
other researchers. For example, Equation (3-41) gives 
ABA which closely fits Hanson's (1971) empirical results, 
and Equation (3-45) gives a ratio of transmittance which 
closely agrees with the empirical results of Budyko (1956) 
and Kimball (1928). Finally, the mean value of ABA in 
clear skies (ABA =.137) corresponds closely to the tropo- 
spheric mean absorptivity of 15 percent reported by London 


C1957 )r 
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IV. TERRESTRIAL RADIATION 


A. THEORETICAL AND EMPIRICAL BASIS 

Sasamori (1968) developed empirical formulas for 
computing emissivities for long-wave flux calculations 
associated with the NCAR General eteeaneeaion Model. 
These empirical emissivity formulas are matched to the 
values developed by Yamamoto (1952). from radiative transfer 
theory, which have also been built into his well-known 
Radiation Chart. This chart has proved to be quite accu- 
rate for computational purposes, and is used here as a 
schematic guide for integration of the radiative transfers 
through the various layers of interest im the soundings 
discussed in Section II. 

For utilization in the FNWC heating package the neces- 
SatvecKiLocivitcy Formulas Of [Sasamorl are those which will 


lead to expressions for 


* 

Flo = net IR flux at earth 
* 

Fe = net IR flux at level k = 6 
* 

F. = net IR flux at level k = 2 


* * 


x * 
and F, - F. will then give the IR 


The differences Pe - Foo - 6 


flux-divergence (1i.e., cooling rates) in the layers (6, 10) 
ana (2, 6)5 respectively, A workable scheme for making 
such computations is presented below for various combinations 


Oimerougecover ([CiEtl), CL(2))]. 
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For example, in order to compute 
25100 
* 
Fio = Bio -f EG (Ue Gc, 1). <d8 (4-1) 
B=0 


one must have a good representation of the emissivity Te 


as a function of both water vapor and CO. absorber masses 


2 
along the sounding, that is along the patheroft integration. 
Figure 7(a) illustrates the concept involved at level 

k = 10, with clear skies, where the hatched area represents 
the downward flux at k = 10. Here the radiative sounding 


is converted into emissivities as indicated in the figure. 


The hatched area denoted by B 
0 


al 
€ dB 
we 
B=0 
is then the key quantity to be evaluated. For this integral 


evaluation, Sasamori (1968) proposes emissivity formula- 


tions of the type 


ean S10) 


- { .240 log,,[U(8, 10)] + 622} 


+ 


Pore G2 fa. BRC COG luce Oj 026617 3) 
[log,, C(8, 10) + 1.064] } aes 


Lemme npetamiie Sad 7210 7k where a 1s very nearly tempera- 
ture independent. ihiewiuanmeley slmethe Lirst brace of Equa= 
tion (4-2) indicates the water vapor emissivity alone, 
whereas that in the second brace is the added contribution 


PyesencCwpdmntmOreche CO. Emissivity transmitted through the 


2 


Weeer Vapor “overlap” near 15 wm. 
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Equation (4-2) was used in the calculation of the atmos- 
pheric emissivities over any limits (6, EO), (4, 10), ... , 
ee lO} £or whieh running mass values U(G@, 10), C(6, 10), 
etc. were available from the radiative sq@undings. Hence 
Equation (4-2) is extremely versatile in its applicability. 
In order to complete the integration from level k = 1 to 
the origin (B = 0, € = 1.0), the U(T,) c(t, ) soundings 
were taken as isothermal at T = Ty over tthe layer (0, 1). 
However, it is still necessary, as Sasamori points out, to 
determine the hatched areas on the Yamamaito diagrams for 
the temperature-dependent emissivity Ee tuo, 10); T,1- 


This was computed after Sasamori as 


e353) 10g, UG@ogel0) =) 244 


Sree LU (0, 10), T,] = 8.34T, 


©-03455=-.705 


<0 10) *18.0/(.353 log,, U(O, 10) + 3.56)] 


+ 07262 { (1. = .625561U(0, 10) + .@@o)°-°) 
*[log,, C(O, 10) + 1.064) } (4-3) 


The final term in the braces of Equation (4-3) involves 
essentially the same form of the CO. corpection to the 
total emissivity as appeared in Equatiom (4-2). The first 
term in Equation (4-3) is, however, temperature-dependent, 
as well as dependent upon the summed absorber masses rela- 
tive to the reference level (taken here as k = 10). It is 


clear that Equation (4-3) is applicable with appropriate 


finaly parameters U(0, k), €(O, k); and Fy. to the 
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determination of the curvilinear portion of the downward 
flux contribution through level k for all of the reference 


level cases to be considered below, e.g., k=2, 6. 


* 


lo’ WITH CLOUDY SKIES 


B. NET FLUX F 
The radiative soundings have been evaluated in the 


form of the parameters U(k, 10), C(k, 10), T level for 


k! 
level. In addition there is listed for each radiative 
sounding at every grid-point the cloud-cover parameters 
CL(1) and CL(2), which in general are not both zero. The 
grid area may then be visualized as comprising the areal 
Praccei1ons W(0, 0), W(l, 1), Wl, 0); and W(0, 1) defined 
by Equation (3-12), of the following cases: 

(a) totally clear skies, (case 0, 0) 

(b) overcast in both layers, eae deel) 

(c) overcast in upper layer only, (case l, QO) 

(d) overcast in lower layer only, (case 0, 1) 

The resultant radiative effect of the true cloud cover 
over the grid area is then reconstructed by multiplying 
the results of the respective cases by the appropriate 
weight factors. It is expedient to discuss first the 
totally clear sky case (Figure 7(a)). 

io) SNet Fr lux F. With Clear Skies 


10 


Tf ene nateched area FY) 6 of Figure 7(a) is evalua- 


ted using the trapezoidal summation rule, then the net 


* 
5 ef 
ux Pio becomes 


oll 








WV pe) eaten A THE) A al ee, Cena s 


e * 
Figure 7(a). Terrestrial net flux Foo with clear skies, 
(case(0,0)). The unhatched area depicts 
Bao 
Fr. 1l-€ d 
10 — Ee crane 
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F = B - -Sfe (8,10) (B, 


10 10 -B,)+(Bo-Be) [e  (6,10) 


6) 


te (8,10) +e (4,10) +e (6,10)] (Be-B,)+(B,-B,) [€ (2,10) 


te (4,10) ] +e, (1,10) te (2,10)1 (BL -B,)+€, £(0,10),7, 1B, } 


B, = 1.170403 * tome oe ly (day)? (4-4) 


Here B, is the Stephen-Boltzman blackbody flux at T 


k k* 


* 
Zee NCt Flux Foo with Overcast Clouds in Both Layers 
Here Figure 7(b) is converted into a numerical 


form uSing the SaSamori type emisSivities [Equations (4-2), 


(4-3)] in the quadrature process. The result is 


& 
Eel, 1) = (8 


10 = B,) ti. a o9€ i C9e 10)] (4-5) 


10 


where By is defined in Equation (4-4). 


* 
Sus Net Flux Foo with Overcast Clouds in Upper Layer 
Only 


Here Figure 7(c) is converted into a numerical 
form by the quadrature process after computing € — using 


Equations (4-2) and (4-3). The result is 


F_..(1, 0) = (B..-B.) {1. - SOLE (Si 10) (B= Be) 


6 


oss 10)+e€, (6, 10)) (B,-B,) 1} (4-6) 


The rationale of Equation (4-6) is that the downward flux 


Moe ccentlallivaB. plus an increment which is in the upper 


6 


ee, 








; * 
Figure 7(b). Terrestrial net flux F 
and low level overcast cloud decks, ( 
unhatched area depicts 
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with both high 
case(1,1)). The 








: * 
Figure 7(c). Terrestrial net flux F 0 with high level 
Overcast clouds in the layer k=4 to Toe only, (case (1,0)). 
The unhatched area depicts 


Bio 


ae. 





right corner of Figure 7(c), given by the summed combina- 


tion or: the products of the form 


as these appear in the brackets in Equation (4-6). 
4. Net Flux Fao with Overcast in Lower Layer Only 
Here the numerical result follows immediately from 
the graphical depiction of Figure 7(d). This case becomes 


a replicate of the result of Section IV(B)2, with 


op §Or1) = (Biot By) fl. ~25€) (9,10) ] (4-7) 


Ds Composite Foo Calculations 
Here we use the principle enunciated in Section IV 
(B), and express Ee ACH pe CL(2)] as the properly cloud=- 
weighted sum of the respective results of cases (a), (b), 


(c), and (dad). The result must have the form [see Equations 


(3-13) and (3-26) ] 


& k 
Foo = W(0,0) Fy 9 60, O) + W(1,1) Fig tit) 


x x 
t+7Wwil, 0) Figo 1,9) WACO, 1.) Fig l9-1) (4-8) 


with thes cloud=weighted fractions, W(0,03, W(1,1), W(1,0) 
and W(0,1) defined in Section III(B)3. The composite re- 


sult simplifies to 


10 10 6 wc 8 


F = (1-cL(2)1 { (B. = BL.) =.5 € (87,10) (By Se en) 
+7 [e, (8,10) ai E (6,10) ) (By = at + 
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Figure 7(d). Terrestrial net flux F 
cast clouds in the layer k=8 to k=9 
The unhatched area depicts 


with low level over- 
tyoevcase (071) ). 


# 
LO 
on 


Bio 
(1 - 


* 
Jie -f eect = 


B=0 


Su 





a 


) 


+ (1=cL(2)) (Len (1) ) | Be7s5 | (c,.(6,10P4e, (4, 10)) (BEB, 


tle, (4,10) +e (2,10) 1 (B,-Bi) +e (2,20) te (1,0) ) (BL -B ) 


sl 


+ E,o(0r10) .7,)] | PemGe) { (B 9-By) [1.-.5e,, (9/1091 


(2-7) 


* 
Note that Equation (4-9) reduces to Fig (O19) when CL(1) = 


* 
CL(2) = 0, and to Fo fil) when CL(1) = CL(2) = 1.0. Simi- 


* 
larly the other two subcases of F may be recovered by 


10 

setting up the weightsS corresponding to CL(1) = 1.0, 

Sbez)) = 0.0 and CL(1) = 0.0, CL(2) = 1.0, respectively. 
This procedure is similarly extemded to levels k=6 


and k=2. The results are as listed in the following sub- 


sections. 


C4 NE FLUX F, WITH CLOUDY SKIES 

Formulation of expressions for the umhatched areas of 
Figures 8(a), (b), (c) and (d), with appropriate weight 
Paceors  W(0O,0);, Will,1), Wii,0), and wW(0 kd eee when 


summed lead at once to the following result: 


x 
Pe =" ({l-CEL(1)) {Bg S LE, (6,8) (BeBe) +E, (4,6) (BE~By) 


Fle (4,6) +e (2,6) )(B,-Bo) +e, (2,6) +e (1,6) ) (B,-B ) 


1 
+ Eye ((0r6) 7) BI Feet her (1) (1—Crigee:) { (By 97Bg) (1- 


5(€,, (6,8) +E, (6,10))1}+ CL (1) {(Bg-B,) (1-5 (e,, (6,8) 1} 


#CL(1) (1-CL(2)) {(B, 9-Bg) (1--5fe,, (6,8) +e, (6,10) 1)} 


(4-10) 
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* 
Figure 8(a). Terrestrial net flux F., with clear skies 
(case (0,0)). The unhatched area depicts 
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* 
Figure 8(b). Terrestrial net flux F., with both high and 
POowelevel overcast layers (case (1,1j)). The unhatched 
area depicts 
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* 
Figure 8(c). Terrestrial net flux F., with high level 
overcast clouds in the layer k=4 to P=6 only (case (1,0)). 
The unhatched area depicts 


ou 








With low level 


* 
Figure 8(d). Terrestrial net flux F_, 
overcast clouds in the layer k=8 to pao Only (case (0,1)). 
The unhatched area depicts 


Equation (4-10) reduces to the properly formulated net 
fluxes for cloud cover cases (0,0), (1,1), (1,9), and (0,1) 


as these are depicted in Figures 8(a),(b),(c) and (da). 


* 
D. NET FLUX PF. WaT CLOUDY SKIES 


* * 
Similar to Fio and F. the net flux at level k=2 can be 


* 
calculated using the trapezoidal rule. The net flux F., 
is schematically depicted by the unhatched areas of Figures 
9(a),(b),(c), and (d). After these areas have been calcu- 


lated anc weighted properly as described in the preceding 


two sections the following formula results 


F-=(l1-cb(1)){B -.5[€ (2,4) (B,-B.)+(€ (2,4)+€ (2,6) ) 
2 8 wc 4 2 wc wc 


(Be-B,)+le. (2,6) +e (2,8) ] (Bo-Be) te, (1,2) (BL-B ) 


1 


#e .((0,2),T:) (8.))]} +(1-CL(1)) Q=cL(2)) 
wc 1 i: 


{(B -B,) [1-.5(€. (2,8)+e (2,10))} 


10 


welt) {By ~-5[€, (2,4) (By-By) +e, (1,2) (B,-B ) 


1 


Fee ((0,2),0,)8, 13 Caen 


This equation also reduces to the net fluxes for the four 
cloud cover cases (0,0), (1,1), (1,0), and (0,1) when the 
appropriate weight factors for the simplified cloud covers 


die "Set ante Equatzon (4-11). 
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F * 
Figure 9(a). Terrestrial net flux F., with clear skies 


(case (0,0)). The unhatched area depicts 


64 








* 
Figure 9(b). Terrestrial net flux F,., with both high and 
low level overcast Cloud decks (case (1,1)). The unhatched 
area depicts 


65 





: 
: 
: 


Figure 9(c). Terrestrial net flux with high level over- 
cast clouds in the layer k=4 to k=6 only (case (1,0)). 
The unhatched area depicts 
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* 
Figure 9(d). Terrestrial net flux F_, with low level 
overcast clouds in the lower layer, f=8 to kK=9, only 
(case (0,1)). The unhatched area depicts 
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E. APPLICATION TO HEAT BALANCE COMPUTATIONS 

In computing the heat balance of the atmosphere, the 
Pavers ©(2;,0) and (6,10) were considered to be undergoing 
solar heating (Section III). Now we require the long- 
wave cooling effects both instantaneously and averaged 
Oven, the day {Section VY). The long-wave cooling in (2,6) 


is symbolized by FU26 which is defined by 
 &* 
FU26 = Fs = /% 


Likewise that in (6,10) is symbolized by FU610 where 


is the flux-divergence in the lower half of the atmosphere. 


For computations of the radiative "balance" at the tropo- 


* 


pause only F 


of Equation (4-11) is required along with 

the net usable insolation. At the surface, a radiative 
balance is computed at each gridpoint involving the absorbed 
insolation and the net flux Fi Of Equation (4-9), the 
latter considered a negative contribution. As will be 

seen in Section V, the long wave flux divergences FU26 and 


FU610 both act as destabilizing influences in their respec- 


tive layers. 


Ee SlLALISLTiICALS RESULTS 
ie Downward Flux with Clear Skies 
Brunt (1932) demonstrated that the nighttime radia- 


tive measurements at the earth were expressible ina 
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statistical form 


Be 3a) le: bve ) (4-12) 


where Bio is the Stephen-Boltzman black body radiation 


4 , 
=o a ) at the surface, e is the water-vapor 


(Bio 10 


pressure in mb, and a and b are the regression coeffi- 


The regression coefficients 


* 
cients, and F, = Bio ~ Fio° 


| 


referred to by Brunt were shown to be subject to wide 


variations depending on Bio and ve. 


In this study a statistical test of Equation (4-12) 
was carried out by linear regression using a program from 
the BMD set of statistical programs (Dixon, 1973). The 


* 
predictand selected for these computations was Foor which 


included the downward contribution of co, as well as that 


of water vapor on the calculation of the net flux at level 


10. The statistical equation that resulted was 


. 
Big ap oc ctOalcies = crete. Ye), R = .995 (A= 13) 


The correlation coefficient of Equation (4-13) shows that 


* 


10 and therefore 


and Ye are Good predictors of F 


= alo Bio 


Of Fae Using the values of a and b in Equation (4-13) and 


Solving fox Py we obtain 


ry 
ll 


q Bio Jama) = w/a 


ny 
ll 


Bi (-69817 + .01314 Ye) (4-745 
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The high value of the correlatiom coefficient indi- 
cates that Equation (4-13) may be regarded as an oceanic 
version of the Brunt F4-equation. Note Rowever that most 
of the dependence of F, is with B and that the variation 


d 10 


* 
of Foo with Biyve is relatively slight. This 1s as to be 
expected over the ocean, where nearly comstant relative 


humidities might well exist over large areas. The mean 


values of Pa and Boo for clear-sky cases are as listed be- 
low 
ee ~ F 
Ps = 637.09 ly day 
Boe] 643.47 ly day a 
10 ; aa 


so that Ee averages 75.5 percent of B These statistical 


10° 
averages are based upon the 57 grid-point soundings in the 


Northern Hemisphere only (those in the Southern Hemisphere 


have been deleted in order to eliminate tropical bias). 


eo potatistical Relationship between Cloudy Sky and 
Clear Sky Cases 


Based upon the radiative sounding parameters listed 
in Table I(b) and Equation (4-9) it was possible to compute 
F, (CL (1),CL(2)] and Fi (0,0) based upon otherwise identi- 
Ga Beene eae 

FOm@ enenpurpose Of this statistical study the ratio 
EG (CE) menny/E 6 (0), 0) was formed as the predictand, 
while an approximation to the composite total opaque cloud 


cover at the grid-point, defined in Equation (3-43) as 


Cen Clg) mete). ( 2) Gl=CrL C1 )) 
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for each cloudy sky case was used as the predictor. The 
simple regression formula based upon a best-fit to the 57 


Northern Hemisphere data soundings became 
4 i, 0 Lo 75 > 
Fig (Ch) = Figo! pO l=. CL) (4-15) 


with a regression coefficient, R= .983. A regression 
analysis using cr? as the Beale tor resulted in only a 
very small improvement of R. The hemispheric mean values 
of the various parameters in Equation (4-15) were as 


follows 


Pig (ch) = ed? 4.8) ly Aces = 
* -1 
Pig (0°0)= 206,398 ly day (4-16) 
CL = 0.6103 


Here the superior bar means the sample mean was taken 
over the 57 Northern Hemisphere grid-points. 

The mean results of (4-16) indicate ee the sur- 
face net flux 1S decreased by the ratio @83.97206.38) = .455 
with CL = 0.6103. This is to say that the downward flux 
al he is on the average 45.5 percent larger than with the 
clear-sky case for the cloud-parameterization of Equation 
(2=13). 


k 


tk 
si F (CL) Relative to F. (0,0) 


Here the statistical parameter sought is the pre- 


da ctand 


* 
ae 

Y = ——— (4-17) 
2 


ripe 





where both numerator and denominator are computed directly 
from Equation (4-11) with (a) the known CL(1), CL(2) used 
in the numerator, and with (b) CL(1) = CL(2) = 0 in the 
denominator. The regression result was cast in the desired 


form 
Y = 1 =- A,(CL) - A (EL) 7 
1 2 


with Ay A, to be determined. The resulting regression 


equation was 
2 
Y= 1 = .0742(CE) = £3450 (CL) (4-18) 


with a multiple correlation coefficient R = 0.9431, anda 
standard error of estimate .0529, The mean values of 


* * 
HEGEL) 7 F, (0, 0), were 
* 
F,(CL) = 453.19 ly day 


* _ 
F,,(0, 0) = 517.70 ly day : 


CL = 0.6803 


—_ 


The mean CL here was based upon the 75 sounding-sample 
which contained a larger proportion of tropical cases than 
did those of the preceding section. 

It is to be noted that the infrapyed loss to space 
has not been greatly impeded by the assumption of normal 
cloud cover since no thick cirrus clouds have been in- 
cluded in the model. With deep cirrus the “greenhouse 


effect" would be increased for gridpoint columns so affected. 
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For comparison with satellite lomg wave flux to 
space, which is usually given in ly afin, the foregoing 
FS values should be divided by 1440 min Pasp oe 

Comparison with satellite measures of long wave 
flux should be made with Fit rather than with F, which 
includes a downward flux contribution from the strato- 
Sphere. For this purpose, Pomc aetons were made with 
Fy of Equation (4-11) but omitting terms involving both 
Revels "2" to "1" and "2" to "Oo", which are always asso- 
Baeed with downward flux (see Figure 9). The comparison 
of the grid-point sounding mean Fi [cL(1),CL(2)] was then 
Made with the NIMBUS III global average long wave flux 


to space for the appropriate date (after Raschke et al, 


LOG ye The comparison is as follows 


F,[CL(1), CL(2)]= .407 ly min” 


long-wave flux to sSpace= .352 ly atin (after Raschke) 


The comparison is quite good considering that the results 
of Raschke et al are Northern Hemisphere means for a two- 
week period in April 1969, whereas the results of this 
study, .407 ly “Sune is essentially based upon a mid- 
Pacitic GCross—section On April 25th, which tends to be 


warmer than the globe as a whole. 
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V. MERIDIONAL CROSS-SECTIONAL DEPICTIGCN OF RESULTS 


In the computer programming of this model many outputs 
have been generated, so, for ease of examination of the 
output, an atmospheric cross~section along each of the 
three meridians studied was developed. This form of pre- 
sentation was designed to give both the imcoming and the 
outgoing radiations at the selected data levels. In addi- 
tion, radiation balances were computed within key atmos- 
pheric layers and at the earth's surface at each gridpoint. 
In order te indicate the relevance of the cloud-amounts 
CL(1),CL(2) and CL to the resulting balances, these values 
were also plotted at appropriate levels (Figure 1l, in 
parentheses) at gridpoints. These cross-sections were 
then interpolated to give the above results at five degree 
latitude gridpoints along the selected meridians. This 
was done for ease in calculation of a composite value of 
each of the computed radiative quantities at identical 
latitudes between 20°S and 65°N, wherever possible. 

Lines 1 and 3 extended from 20°S to 55°N and had identical 
latitude grid-points. Line 2 extended from the equator to 
65°N and therefore permitted extension of the cross-section 
Pemos lass Nclatzttudges South of the equator only two values 
were available for the composite computations and at lati- 


tudes higher than 55°N, only one value was available. 
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Between these extremes all three meridional cross-sections 
were used in the calculation of the composite case, Fig- 
ure ll(a) and (b) depicts the composite cross-section with 
Figure 10 eel tite key for locating the desired values. 
In the cross-section to be depicted the insolation 
rates at level k=2, (QAVG) are given as averages over the 
24-hour day (25 April 1973). This was done by employing 


Equation (3-1), but replacing Cos Z by Cos Z where 
Gos 2 = [h Sind Sind + Cosd@ Cosé Sin h]/T ~ (5-1) 


and h 


hour angle = Arc Cos(-Tand Tanéd) 


latitude 


oOo & 
I! I 


declination (Smithsonian Tables (List, 1963) ) 


Here the double bar over Cos Z represents a time mean over 
a 24 hour period, that is for a time period compatible 
with the long-wave flux transfers. 

At each latitude a value of the instantaneous absorp- 
tivities, reflectivities or transmissivities was known for 
the time OOOOGMT, and for the purpose of the cross-section 
presented here, these were assumed constant even though 
the insolation rates changed over the day. The IR trans- 
fer rates and their flux divergences were assumed to be 
constant both for the day and the composite grid-point as 
determined from their OOOOGMT soundings. Thus it was 


possible to compute: 
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a) the earth-atmosphere radiative balance at level 


k=2 by the formula 
, te 
QTROP = QAVG - (REFF + F,) (5-2) 


b) the radiative balance in the layer (2, 6) as given 
by 


NET26 = ABS26 + FU26 (5-3) 


c) the radiative balance in the layer (6, 10) as given 
by 


NET610 = ABS610 + FUG61LO (5-4) 


d) the radiative balance at the earth's surface as 
given by 


* 
NETSFC = ABSFC + Pio (5-5) 


% 


In addition to these mean grid=-point items the cloud covers 
CL(1), CL(2) and (CL) in parentheses are depicted in the 
applicable layer after averaging to the grid-point. At the 
tropopause, and at the earth's surface a global average of 
radiative balance or net was obtained by cosine weighting 
as follows 


a >, Q(o.) Cos @. 


yy Cos >. 


where Q(¢,) is the radiative net value under consideration 


ated ac i Cude >.- Some of these values are contrasted with 
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the other climatological results cited below 


This Study Other Data 


1) global balance at 
the tropopause - 048 ~O61, Raschke 


2) global balance at 
the surface 2166 -140, Budyko 

The outgoing long wave radiation in these computations was 
also in excellent agreement with that of Raschke et al 
(1973) as noted in Section IV. The global balance of total 
radiation at the tropopause is less than that reported by 
Raschke et al for the same general time of year. However, 
it was noted that this result primarily reflects the 
smaller value of planetary albedo (ALB = .300) reported 
from satellite observations, 

Another set of results for comparison with those calcu- 
lated here are due to Budyko (1956). Budyko calculated 
the radiation balance at the surface for ten-degree lati- 
tude bands. In his study, these values were taken for the 
balance at mid-points of these bands and compared with the 
computed output from the present model. The following 


table shows this comparison. 


Budyko 
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All of the above values are in cal cm min a These 

comparisons would have been in closer agreement if the 
results of this study had been in the form of an annual 


average as Budyko's results were, 
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Vi. CONCLUSIONS 


Using three meridional lines of soundings over the 
Pacific Ocean a radiative model for inclusion in the 
FNWC numerical weather prediction scheme was evaluated. 
These three meridians of data-soundings were used to 
compute a space-time average of the radia@tional para- 
meters for one day over the Pacific Ocean. This space- 
time average agreed fairly well with values derived from 
other studies, including satellite measurements. There 
were some Minor differences which could be ascribed to 
the fact that the data used in this experiment was a one 
day space-time average while the comparison data was an 
average of several weeks for the satellite data, and 
over years in the case of surface climatological data, 

Differences could also be induced by the parameteri- 
zation of cloud amounts as well as by the treatment of 
cloud reflectivity and absorptivity for solar wavelengths. 
In this model, only two layers of clouds were considered 
and the same empirical values of the reflectivities and 
Deere were used for all latitudes, which could 
cause discrepancies when compared with satellite data. 

Beerecen the results of this data treatment agree 
Paacly  weliawten climatological Comparisons cited at the 
end of Section V, refinements may be possible in the 


present model. GBhis. model should be Cested using data 
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at different seasons of the year and over varying 
locales. With a full year of sampling, refinements in 
the cloud parameterization, atmospheric absorption 
values and reflection treatments could be developed 


for accurate use in the FNWC predictive model. 
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